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j^. . We study the properties of a future singularity encountered by a 

CN \ perfect fluid observer in tilting spatially homogeneous Bianchi cos- 

mologies. We derive the boost formulae for the Weyl tensor to estab- 
lish that, for two observers that are asymptotically null with respect 
to each other, their respective Weyl parameters generally both tend 
to zero, constant, or infinity together. We examine three classes of 
00 ■ typical examples and one exceptional class. Given the behaviour of 

| the Weyl parameter, we can predict that the singularity encountered 

is a Weyl singularity or a kinematic singularity. The analysis suggests 
that the kinematic variables are also useful in indicating a singularity 



q-i in these models. 



PACS numbers: 98.80.Jk, 04.20.-q 



1 Introduction 

A cosmological model (Ad,g,u) is defined by specifying the spacetime ge- 
ometry, determined by a Lorentzian metric g defined on the manifold Ai, 
and a family of fundamental observers, whose congruencd^ of world-lines is 
represented by the 4-velocity field u. The covariant derivative u a - b of the 
4-velocity field is decomposed into kinematic variables according to 

U a -b = &ab + ^ab + H(g ab + U a U b ) - U a U b , (1) 

where a ab is the rate of shear tensor, uo ab is the vorticity tensor, H is the 
Hubble scalar, and u a is the acceleration vector. Sometimes there is another 



1 A congruence is a family of curves such that through each point there passes precisely 
one curve in this family. 
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preferred family of fundamental observers, u, whose covariant derivative de- 
fines a corresponding set of kinematic variables (a a b,uj a b, H,u a ), which are 
different from the first set. 

The study of spatially homogeneous (SH) Bianchi cosmologies forms part 
of the effort to understand the underlying dynamics of anisotropic and inho- 
mogeneous universes. In recent papers [U [2] we examined the future asymp- 
totic dynamics as experienced by the perfect fluid observer in Bianchi cos- 
mologies with a tilted perfect fluid with linear equation of state p = (7 — 
focussing on physical properties of the models as seen by the fluid observer. 
We found that given a Bianchi type, for a large enough value of the equation 
of state parameter, 7, the perfect fluid observer can encounter a singularity, in 
spite of the fact that the SH observers, moving geodesically and orthogonally 
to the spatially homogeneous hypersurfaces, do not encounter any singularity 
into the future [3]. Although there is no universally accepted definition of a 
singularity [4, Section 9.1], we shall make use of the congruence in defining 
a singularity within the context of a cosmological model: a singularity is en- 
countered by a congruence if the proper time is finite and inextendible. This 
definition is consistent with the general definition adopted in the early work 
on Bianchi cosmologies by Ellis & King, which states that a singularity is 
recognized by the existence of an inextendible curve which has a finite length 
as measured by a generalized affine parameter [5j page 1211. E The encounter 
of a singularity is accompanied by an extreme tilt limit (i.e., the fluid ob- 
server become asymptotically null with respect to the SH observer). The 
singularity encountered by the fluid observer has the following properties]! 

• The proper time needed to reach the singularity is finite, 

• The Hubble scalar H and some other kinematic quantities diverge, 

• The matter density tends to zero (i.e., the singularity is not a matter 
singularity) 



2 Many also use the stricter definition of a singularity as geodesic incompleteness. How- 
ever, in our case this is not satisfactory since the fluid does not follow geodesies; thus we 
will not require the congruence to be geodesic. 

3 Here, we refer to the properties of variables associated with the orthonormal frame of 
the fluid observer. 

4 Since the source is a perfect fluid with a linear equation of state, this means that all 
components of the Ricci tensor also tend to zero. 
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All components of the Weyl curvature tensor converge for some cases and 
at least one component diverges to infinity for other cases. In this paper we 
shall further examine the behaviour of the components of the Weyl tensor 
(according to the fluid observer), with the aim to predict whether they con- 
verge or diverge. In order to do so, we recall the following concept regarding 
the Weyl curvature tensor. 

The Weyl parameter 

The Weyl tensor is defined as the trace-free part of the Riemann curvature 
tensor and, consequently, is not directly involved in the Einstein field equa- 
tions. In some sense, the Weyl tensor describes the 'free gravitational field' 
and therefore is of particular interest. Penrose [6] suggested that the Weyl 
tensor is related to a measure of a 'gravitational entropy' and therefore could 
be used to shed light on the initial state of the Universe. In this context, 
the Weyl tensor and its Weyl scalars have been used to characterise different 
behaviours of non-tilted Bianchi models pJIEllS]- 

Wainwright, Hancock and Uggla [TU] studied the late-time dynamics of 
Bianchi type VIIo cosmologies with a non-tilted perfect fluid. They intro- 
duced a quantity W called the 'Weyl parameter' (see [10J, p 2580): 

W = -^, where W 2 = \{E ab E ah + H ab H ab ). (2) 

They coined the term0 Weyl curvature dominance' to describe the phe- 
nomenon in which W — > oo. This means that the ratio of some compo- 
nents of the Weyl curvature tensor and the square of the Hubble scalar tends 
to infmityl§ This phenomenon also occurs in Bianchi type VIIo and VIII 
cosmologies with a tilted perfect fluid [HI [T2] . 

In [TOj [TTJ [12], the Weyl parameter is defined with respect to the SH 
observer. However, the Weyl parameter can also be defined with respect 
to different observers, and particularly the fluid observer in tilted Bianchi 
cosmologies 

5 In [9] the term 'extreme Weyl dominance' was used for the case W — > oo while 'Weyl 
dominance' was used for the case when the Weyl curvature invariant dominated the Ricci 
invariant but W was bounded. 

6 It is worth noting that both W and H tend to zero in [TU] , 

7 The Weyl parameter, unlike the invariant Weyl scalars, is an observer-dependent quan- 
tity. 
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Although the limits of the Weyl parameter along two congruences of 
worldlines are in principle different, we conjecture that in general they both 
tend to zero, constant, or infinity together. We shall refer to this as 'the 
Weyl parameters have the same convergence/divergence property'. The ar- 
gument is considered in two cases: whether or not the two congruences are 
asymptotically null with respect to each other (v 2 — > 1). 

If the two congruences are not asymptotically null with respect to each 
other (v 2 -f* 1), then T = (1 — v 2 )~ 1 ^ 2 is bounded. Then, from the boost 
formulae (134|) - (!35|) for the Weyl components, in general the Weyl compo- 
nents Cabcd and Cabcd relative to two different congruences have the same 
convergence / divergence property. Also recall from the boost formula for the 
Hubble scalar H (pQ, equation (B.15), reproduced below) 
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e^IV) - JL. v ^(r) + 3TH - 2Ya^ + Vii^ 



(3) 



that H and H also have the same convergence/divergence property. It then 
follows from (T5]) that generally the Weyl parameters W and W have the same 
convergence/divergence property. 

In the more interesting case for this analysis, when the two congruences 
are asymptotically null with respect to each other {v 2 — » 1), then r — ► oo 
and, in general, C^fccd becomes of order T 2 Cabcd- From we see that H is 
of order TH as T — > oo. As a result, generally we have 

Cabcd . r i Cabcd / A \ 

— — is of order — ; (4) 

H 2 H 2 

i.e., the Weyl parameter VV of a second congruence is generally of the same 
order as the Weyl parameter W of the first congruence, and they therefore 
have the same convergence/divergence property. Exceptions are possible 
when cancellation occurs in the leading order terms in all the Weyl compo- 
nents, as we will show occurs in the LRS Bianchi type V cosmologies. Indeed, 
we will show that, when the argument works, we can use it to predict the 
type of singularity encountered by the fluid observer. 



2 Terminology 

We first review some terminologies about singularities and the asymptotic 
dynamics of the Weyl tensor. 
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There are several classifications of singularities in the literature. We are 
concerned with one particular classification regarding the convergence/diver- 
gence of the Weyl curvature tensor at the singularity. 

Weyl singularity 

Collins and Ellis [131 p 88] used the terminology 'Weyl singularity' to de- 
scribe the following: at least one component of the Weyl curvature tensor 
(with respect to the orthonormal frame of the observer who encounters the 
singularity) diverges]^ 

On the other hand, the matter density and all components of the Weyl 
tensor can converge as the singularity is approached, while some kinematic 
quantities (most importantly, the Hubble scalar H) diverge. We shall call 
such a singularity a 'kinematic singularity' when it is not a Weyl singularity 
or a matter singularity. In other words, relative to a congruence, a kine- 
matic singularity is characterized by the blow-up of one or more kinematic 
variables in finite proper time, while all components of the Weyl and Ricci 
tensors remain bounded. This terminology is new. A prototypical exam- 
ple of a kinematic singularity is the initial singularity encountered by the 
fundamental observers of the Milne universe p3] . 

Later, in section HI we shall discuss this classification in the literature. 

Weyl blow-up 

The concept of Weyl singularity combines two phenomena - the blow-up of 
Weyl and the occurrence of a singularity. To separate them, we introduce the 
terminology 'Weyl blow-up' to describe the phenomenon when some compo- 
nents of the Weyl tensor diverge asymptotically with respect to a particular 
observer, regardless of the occurrence of singularity]^ Weyl blow-up can be 
simply stated as 

M/^oo, (5) 

where W is defined in The phenomena W — *> const, or W — > will be 
called 'Weyl convergence'. This terminology is new. We have not encountered 

8 Collins and Ellis actually used the terminology 'conformal singularity', but we feel 
'Weyl singularity' is more appropriate and avoids confusion with the 'conformal singularity' 
in the context of isotropic singularities [TH [15] . 

9 It is possible for Weyl blow-up to take an infinite proper time to occur, in which case 
there is no singularity. 
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examples in tilted SH cosmologies in which W is bounded but whose limit 
does not exist. 

The Weyl scalars may or may not converge asymptotically when Weyl 
blow-up occurs for one observer. There are three scenarios. In the trivial 
case where the Weyl tensor is identically zero, we have Weyl convergence 
for all observers. In the second case, if at least one of the four Weyl scalars 
diverges, then we have Weyl blow-up for all observers. In the third case, if all 
four Weyl scalars converge, then one observer may experience Weyl blow-up 
while another does not. The third case is the case of interest here. 

Future singularity and Weyl blow-up 

In the next section we shall examine tilted Bianchi cosmologies that have 
an extreme-tilt sink. On approach to the sink, the fluid observer may or 
may not encounter a future singularity, and may or may not experience Weyl 
blow-up. We summarize the four possible scenarios with a Venn diagram 
(see figure [Q . 



/ Weyl 1 


?low-up \ 






/ No 


Weyl 


Kinematic 




1 singularity 


singularity 


singularity 






\ Singularity / 


No singularity or Weyl blow-up 



Figure 1: Four possible scenarios for the fluid observer at late times. 

We shall use the Weyl parameter to heuristically predict whether Weyl 
blow-up occurs if the fluid observer encounters a future singularity; i.e., 
whether the singularity is a Weyl singularity. We will also give a special 
example in which the prediction breaks down. 
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3 Typical Examples 



In this section we first examine three typical examples of tilted Bianchi cos- 
mologies that have an extreme-tilt asymptotic regime at late times. For more 
details of the asymptotic dynamics of the tilt variable, the fluid cosmological 
time, and the fluid Hubble scalar, see [U Sections 2.1, 3.2.1, 3.2.2 and 5]. 
The three examples illustrate different limits of the Weyl parameter of the 
fluid observer: 



W -> 0, W -> oo, and W -> const. 



respectively. 

In all of the examples we have used the time parameter r to be the 
dynamical time of the SH observer; i.e., r is defined as 

where t and H is the cosmological time and the Hubble parameter of the 
SH observer, respectively. We shall be using r in the asymptotic expansions 
of both SH and fluid quantities. For those who would prefer to express 
the expansions in terms of the fluid cosmological time £fl U id, simply use the 
relation (valid along each fluid worldline but not across the fluid congruence) 



"tfluid — 7} O.T. 

tl 

For further details and discussion see (TJ page 3575]. In pQ we claimed that 
the limits of the kinematic variables along the SH congruence (or the normal 
congruence in the inhomogeneous example) and those along the fluid con- 
gruence are the same. We take this opportunity to justify that claim. The 
examples in [TJ are special in that they have spatially homogeneous limits 
along the SH or normal congruence as r tends to infinity. The fluid con- 
gruence is also special in that it must pass by all events passed by the SH 
or normal congruence. Thus, moving along a fluid worldline, the limit of 
a kinematic variable must be the same as the limit along an SH or normal 
worldline. 

Irrotational Bianchi type V 

The dynamics of irrotational tilted Bianchi type V cosmologies was studied 
by Hewitt and Wainwright [T7]. As the solutions approach the extreme-tilt 
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Milne solution M (which is a sink for 7 satisfying | < 7 < 2) at late times, 
the (orthogonal) SH observer sees 

H^H e- T , W « # 2 |£_o|e- 4T ', W « |£- |e~ 2T , (6) 

as r — ► 00 As a result 

#^0, W^O, W^O, (7) 

for all § < 7 < 2. 

According to the (tilted) fluid observer, however, the asymptotic expan- 
sions are (from the boost formulae (j3J) and flMl) - fl35l) ) 

r-r.-p(^). (8 ) 

W - 2v/2r^ 2 |S_ | exp ( 2(7 2 7 _~ 7 10) ^ , W - ^(2 - 7) 2 |£-o|e^, 

(9) 

as r — ► 00, so 

# -> 00 for § < 7 < 2, 

# -> 00 for f < 7 < 2, (10) 
W -f for all § < 7 < 2. 

We notice that the convergence/divergence property of H and W is observer- 
dependent, while W tends to zero for both observers. For the fluid observer, 
the 7 threshold for encountering a singularity is lower than that for Weyl 
blow-up. A singularity is not necessarily a Weyl singularity in this case (i.e., 
it could be a kinematic singularity). 



Bianchi type VIIo 

The dynamics of tilted Bianchi type VIIo cosmologies was studied by Hervik 
et al [UJ. As the solutions approach the extreme-tilt limit P4 (which is a 

10 The asymptotic expansions for kinematic and matter variables are obtained by lin- 
earizing equations (2.10) and (2.11) in [17], and solving them. These are then substituted 
into equations (14) and (15) in [T5J Appendix 3] to obtain the asymptotic expansions for 
the Weyl components. 
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future attractor for 7 satisfying | < 7 < 2) at late times, the SH observer 



4 
3 

sees 

-2t ti/ _ . „— 3r 



H~e~ 2T , W~ e - 3r , W~e T , (11) 
as r — ► 00. As a result, 

#^0, W^O, W^oo, (12) 

for all § < 7 < 2. 

According to the fluid observer, however, the asymptotic rates are 

r ~ eX p (yr^) > # ~ ex P (yr7 r ) ' (13) 

# ~ exp ( 14 r ) , W ~ e T , (14) 



as r — > 00, so 
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— > 00 for § < 7 < 2, 



W -> 00 for f < 7 < 2, (15) 
W -> 00 for all § < 7 < 2. 

Again, we notice that the convergence/divergence property of i7 and W is 
observer-dependent, while W tends to infinity (i.e., Weyl dominance) for 
both observers. For the fluid observer, the 7 threshold for encountering a 
singularity is higher than that for Weyl blow-up. A singularity is necessarily 
a Weyl singularity in this case0 



Bianchi type VII^ 

The dynamics of tilted Bianchi type VII^ cosmologies was studied in [201 12] • 
As the solutions approach the extreme-tilt vacuum plane wave £_(VII/ l ) 
(which is a sink for 7 satisfying 5+ ^ E < 7 < 2, — | < E + < 0) at late times, 
the SH observer sees 

H ss # e (1_2S+)r , W ~ # 2 , (16) 

as r — > 00. As a result 

H^O, W^O, const., (17) 

11 Tilted Bianchi type VIII cosmologies also exhibit the same qualitative behaviour (see 

mm)- 
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for a11 s+fc < 7 < 2. 

According to the fluid observer, however, the asymptotic rates are 



H 



~ exp 



r 




H 2 



(18) 



(19) 



as r 



oo, so 



H 

W 

vv 



oo 



oo 



const 



for |(2-S + ) < 7 <2, 
for |(2-S + ) < 7 <2, 
for all 5^ < 7 < 2. 



(20) 



Again, we notice that the convergence/divergence property of H and W is 
observer-dependent, while W tends to a constant for both observers. For the 
fluid observer, the 7 threshold for encountering a singularity is equal to that 
for Weyl blow-up. A singularity is necessarily a Weyl singularity in this case. 

From the three typical examples above, we have formulated the following 
conjecture. 

Conjecture: In general (except in a zero-measure set of special solutions 
where the leading order term of each Weyl component vanishes) , the conver- 
gence/ divergence property of the Weyl parameter W is observer-independent. 
In the extreme-tilt limit, according to the fluid observer, 

• If VV — > 0, then the 7 threshold for encountering a singularity is lower 
than that for Weyl blow-up. 

• If VV — > 00, then the 7 threshold for encountering a singularity is higher 
than that for Weyl blow-up. 

• If VV — > const., then the 7 threshold for encountering a singularity is 
equal to that for Weyl blow-up. 

This conjecture works for typical examples, but may fail due to simplifi- 
cation in special solutions. Let us present one such example. 



10 



Exceptional example: LRS Bianchi type V 

The LRS Bianchi type V cosmologies are a special case of the irrotational 
tilted Bianchi type V cosmologies. They were studied by Collins and El- 
lis PJ. 

The shear component £_ is zero in this class, and the SH observer does 
not see the asymptotic rates (Q, but rather 

H ~ H e~ T , W^^H 2 \Z +0 \e- 6T , W « ||S +0 |e- 2 ^, (21) 

as r — » oo. Nonetheless, (J7J) still holds: 

H^O, W^O, W^O, (22) 

for all | < 7 < 2. 

It turns out that there is only one independent component of the Weyl 
tensor C a bcd, namely 

E + = {H + a + )a + + \^-V 2 . (23) 

The boost formulae (I34l) - f[3"5"j) then imply that E + is the only nonzero com- 
ponent in any other frame, and that 

E + = E +] (24) 

i.e., the leading order terms T 2 C a b c d in Cabcd cancel. As a result, equation (J4]) 
fails, and this is a special example when general behaviour predicted by the 
conjecture fails. 

Indeed, the fluid observer sees 

r^Toexp^^^V) , (25) 



W = -# 2 |£+o|e- 6r , (27) 
W - |(2 - 7 ) 2 r 2 |S +0 | exp ( ~ 2 fl~ 2) r ) > ( 28 ) 
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as r — > oo. As a result, 

H — > oo for | < 7 < 2, 

# -> for all | < 7 < 2, (29) 

W -> for all | < 7 < 2. 

Comparing with equation (jlOp . the quantity W never tends to infinity. 

4 Discussion 

The future asymptotic dynamics in tilting SH Bianchi cosmologies has re- 
cently been studied, with an emphasis on the physical properties of the mod- 
els as experienced by the fluid observer (Tj [2]. In this paper we have been 
primarily interested in the properties of the singularity that the perfect fluid 
observer can encounter, with particular emphasis on the behaviour of the 
Weyl tensor (according to the fluid observer). As we have seen in the exam- 
ples, Weyl blow-up is not a good indicator of a singularity. A better measure 
is the blow-up of kinematic variables. 

We now discuss the work by Collins and Ellis p3] and Ellis and King [5] 
on Weyl blow-up at the singularity. In Collins and Ellis [13], it is claimed 
that the future singularity in the LRS Bianchi type V example is a Weyl 
singularity (i.e., a 'conformal singularity', see page 97), in that some compo- 
nents of Weyl tensor diverge. This is not the case, for the only component of 
the Weyl tensor is E + , and it tends to zero. All components of the Riemann 
tensor tend to zero at the singularity. Some kinematic variables diverge, and 
are responsible for the singular behaviour. 

Ellis and King [5] classify singularities as 'curvature singularities', 'locally 
extendible singularities' or 'intermediate singularities'. Curvature singulari- 
ties correspond to our Weyl singularities with the blow-up of at least one of 
the Weyl scalars; locally extendible singularities correspond to our kinematic 
singularities; and intermediate singularities correspond to our Weyl singular- 
ities with bounded Weyl scalars. According to this classification, all of our 
examples of Weyl singularities are intermediate singularities. We prefer the 
terminology 'kinematic singularity' over 'locally extendible singularity' since 
in our examples the fluid worldlines are certainly not extendible at the sin- 
gularity, due to the blow up of the kinematic variables; indeed, our examples 
are not 'locally extendible' in the sense of Clarke [22] . 
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Conclusion 

In discussions regarding singularities, usually attention is focussed on the 
Weyl tensor components and the matter density. We have argued that in a 
cosmological context more emphasis should be put on the kinematic variables. 

In this paper we have derived the boost formulae for the components of 
the Weyl tensor. We have applied the boost formulae to establish that, for 
two observers that are asymptotically null with respect to each other, their 
respective Weyl parameter generally both tend to zero, constant, or infinity 
together. We examined three classes of typical examples and one exceptional 
class. Given the behaviour of the Weyl parameter, we can predict that the 
singularity encountered is a Weyl singularity or a kinematic singularity. The 
examples suggest that the kinematic variables are more useful than the Weyl 
tensor components in indicating a singularity. 
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Appendix: Boost formulae for the Weyl tensor 
components 

In pQ the boost formulae for the kinematic and the matter variables are 
derived. Here we extend the formulae to include the Weyl components. 

Recall from [TJ that the boost formulae for the orthonormal frame vector 
fields are: 



e = Te + Tv^ 
e a = Tv a e + B^e^ 



where 



5^ + 



r + i 



-v n v 



1 



(30) 
(31) 

(32) 



The Weyl curvature tensor C a b c d can be decomposed into 'electric' and 'mag- 
netic' components as follows: 



a 



aO/30 



E, 



a/3 , 



fJ,V > 



/3 7 



7M ■ 



(33) 
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Boosting the orthonormal frame results in the following boost formulae for 
the Weyl components: 

2r 2 f2r + 1) r 4 

E aP = (2r 2 - l)E af3 ^ + 1 V ^, + -^—^EprWv^vn 

2r 3 

+ 2T 2 v^\ a H p)v - ^—v^-^v^y (34) 

2T 2 (2T + \) T 4 
H aP = (2r 2 - l)H aP L-JL-LvPH^vq + -^—^H, u v^v {a v p) 

2r 3 

- 2V 2 v^\ a E p)v + ^—v^^E^y. (35) 
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